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An evaluation of upper troposphere NO, with two models 

Joyce E. Penner , 1 Daniel J. Bergmann , 2 John J. Walton , 1 Douglas Kinnison , 2 Michael 
J. Prather , 3 4 5 Douglas Rotman,^ Colin Priced Kenneth E. Pickering, and Steven L. 
Baughcum 6 

Abstract. Upper tropospheric NO, controls, in part, the distribution of ozone in this greenhouse- 
sensitive region of the atmosphere. Many factors control NO, in this region. As a result it is 
difficult to assess uncertainties in anthropogenic perturbations to NO from aircraft, for example, 
without understanding the role of the other major NO, sources in the upper troposphere. These 
include in situ sources (lightning, aircraft), convection from the surface (biomass burning, fossil 
fuels, soils), stratospheric intrusions, and photochemical recycling from HN0 3 . This work 
examines the separate contribution to upper tropospheric "primary" NO, from each source category 
and uses two different chemical transport models (CTMs) to represent a range of possible 
atmospheric transport. Because aircraft emissions are tied to particular pressure altitudes, it is 
important to understand whether those emissions are placed in the model stratosphere or 
troposphere and to assess whether the models can adequately differentiate stratospheric air from 
tropospheric air. We examine these issues by defining a point-by-point tracer tropopause in order 
to differentiate stratosphere from troposphere in terms of NO, perturbations. Both models predict 
similar zonal average peak enhancements of primary NO, due to aircraft (=10-20 parts per trillion 
by volume (pptv) in both January and July); however, the placement of this peak is primarily in a 
region of large stratospheric influence in one model and centered near the level evaluated as the 
tracer tropopause in the second. Below the tracer tropopause, both models show negligible NO, 
derived directly from the stratospheric source. Also, they predict a typically low background of 1- 
20 pptv NO, when tropospheric HNO 3 is constrained to be 100 pptv of HNO 3 . The two models 
calculate large differences in the total background NO, (defined as the source of NO, from 
lightning + stratosphere + surface + HNO 3 ) when using identical loss frequencies for NO,. This 
difference is primarily due to differing treatments of vertical transport. An improved diagnosis of 
this transport that is relevant to NO, requires either measurements of a surface-based tracer with a 
substantially shorter lifetime than or diagnosis and mapping of tracer correlations with 

different source signatures. Because of differences in transport by the two models we cannot 
constrain the source of NO, from lightning through comparison of average model concentrations 
with observations of NO*. 


1. Introduction 

Tropospheric NO* (= NO + N0 2 ) plays an important role in 
determining tropospheric ozone concentrations. Because 
ozone in the upper troposphere and lower stratosphere is 
especially effective as a greenhouse gas [Fishman et ai, 1979; 
f nr is et al. t 1990], it is important to assess whether 
concentrations of NO* in this region may be altered by human 
activity. Models have demonstrated that the net ozone 
production rate in this region may increase by roughly a factor 
of 2 for NO* increasing from 50 to 200 parts per trillion by 
volume (pptv) [Friedi et ai, 1997]. 
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NO* emissions from aircraft are of particular concern for 
perturbing this region because their emissions are deposited 
directly into the upper troposphere and lower stratosphere 
[Thompson et al ., 1996; Friedi et al . , 1997]. However, the 
background concentration of NO* in the upper troposphere i s 
controlled by a variety of other sources; (1) NO* produced in 
lightning flashes, (2) NO* produced in the stratosphere from 
oxidation of N 2 0, and (3) surface-based sources (soil microbial 
activity, fossil fuel burning, and biomass burning). 

Unfortunately, the source strength for lightning NO*, in 
particular, is not well constrained. O: * might hope to 
constrain it by comparison of predicted NO* with 
measurements. Unfortunately, this is not possible because the 
lifetime for NO* is short and the available data are sparse 
t Emmons et al . , 1997]. Though comparison with existing data 
provides some measure of verification, ideally, one would like 
to also find unique methods to validate the concentration of 
NO* from specific sources because the observed total NO* can 
be derived from many different combinations of sources. Here 
we design a model experiment to isolate and hence identify the 
contribution of individual sources of NO* in the upper 
troposphere as well as the role of transport as represented in 
our two models. 

Transport processes affecting the concentrations of NO* i n 
the upper troposphere/lower stratosphere are represented in 
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models with considerable uncertainty, in part due to failure to 
resolve the relevant spatial scales. The highest-resolution 
chemical transport models in use today are typically T42 
(approximately 2.5* X 2.5°), but even models with this 
resolution cannot resolve the horizontal scales of convection 
that bring NO* to the upper troposphere from the surface. 
Further, they have difficulty representing tropopause fold 
events [Rood et al ., 1992; Holton et aL , 1995] that bring 
stratospheric NO* and O 3 into the troposphere. Here we use two 
tracer transport models as examples to define a range in upper 
tropospheric/lower stratospheric NO* that results from 
differences in treating transport in this region. A similar 
model intercomparison study was recently undertaken by Van 
Velthoven et al. [1997]. Here our results differ, because we 
consider the role of each source in determining the NO* 
concentration and because we use a more realistic treatment of 
the lifetime for NO*. 

We define our experiment in terms of the major classes of 
'primary NO*," for example, NO* that has not undergone 
transformation to the longer-lived species HNO 3 (or 
peroxyacetyl nitrate (PAN)). Source strengths associated with 
surface-based sources (fossil fuel burning, microbial activity 
in soils, and biomass burning) appear to be reasonably well 
defined (Table 1; see also Lee et al. [1997]). Likewise the 
source of total odd-nitrogen (NO v = NO* + HN0 3 + others, see 
below) from the reaction of N 2 0 with O^D) in the stratosphere 
is well defined and thus limits the stratospheric source of NO*. 
Current ’’best guesses" for the lightning source of NO*, 
however, differ by more than a factor of 5 with a total NO* 
possible range estimated as 1 to 25 Tg N/yr [Lawrence et al 
1995; Levy et al ., 1996; Price et al. y 1997a, b]. We also 
examine the contribution to NO* from "secondary NO* " 
derived from in situ photochemical recycling of HNO 3 . 
Concentrations of HNO 3 in the upper troposphere are also 
poorly constrained by measurement, and current models do a 
poor job of representing the available data in any case [Friedl 
et al. y 1997]. Below we show that NO* produced from HNO 3 in 
the troposphere is small over the range of constrained HNO 3 
concentrations from 50 to 200 pptv. 

The upper troposphere and lower stratosphere are obviously 
different regimes in terms of the sources and buildup of NO*, 
and to evaluate the effects of aircraft emissions, they must be 
clearly separated in the model. The chemical lifetime for NO* 
in the upper troposphere is 5 to 10 days, comparable to the 
timescales for venting of boundary layer NO* via convection 
and turnover by large-scale air flow. Thus a balance of 
transport and chemistry determines the distribution of NO* 
from different sources here. In the lower stratosphere, 
timescales for transport are much longer, and a source 
introduced into this region can build up to concentrations 
much higher than those in the upper troposphere. Given the 


longer tr insport times in the lower stratosphere, the chemistry 
will approach a steady state between NO* and the other forms 
of NOy. To evaluate the effects of different sources in the upper 
troposphere and lower stratosphere, it is critical for each 
model to define the surface below which air is essentially 
tropospheric and turns over on timescales less than the NO* 
lifetime. We term this surface the model's "tracer tropopause" 
and define it at each time step and location on the basis of 
simulatiois of downward moving tracers. This surface does not 
necessari y correspond to the model s tropopause defined from 
the temperature lapse rate, but it does accurately separate those 
regions in the chemical transport model (CTM) where sources 
may builc up and NO* should be in photochemical steady state 
from the mixed, non-steady state, multisource environment of 
the trope sphere. The use of this diagnostic allows us to 
evaluate where emissions specified by pressure altitude, such 
as aircraf emissions, are introduced in relation to the modeled 
"tracer topopause." 

In the following we first define our approach, including our 
simplifiec treatment of NO* chemistry, the sources specified in 
the simulations, and a short description of the models. Section 
3 discusses our derivation of the "tracer tropopause" in each 
model and presents results for the concentration of NO* derived 
from the tratospheric source of odd nitrogen. Section 4 shows 
results derived for the predicted NO* concentrations from the 
different tropospheric sources, including the source from 
recycling of HNO 3 . Finally, section 5 presents our 
conclusions, defining for this pair of CTMs the range of 
possible upper tropospheric NO* concentrations associated 
with transport, lightning, and HNO 3 recycling. 

2. Model Description 

The experiments described here used the GRANTOUR/ 
Community Climate Model 1 (CCM1) model [ Penner et a/., 
1991] and the IMPACT model. The latter model is under 
developm< nt as a code capable of using massively parallel 
computer architectures, and the simulations reported here were 
run on 64 processors of the Cray T3d. For this study, IMPACT 
used meteorological fields from the Goddard Data Assimilation 
Office (DAO) stratospheric model that were interpolated to 25 
vertical layers covering the year 1992 [e.g., Schubert et al ., 
1993] (for another paper using these meteorology fields, see 
Weaver et al [1996]). The vertical grid is a hybrid coordinate 
system nude up of 11 terrain-following sigma levels (mean 
pressure levels of 921, 771, 648, 547, 447, 346, 273, 223, 
187, 162, ind 139 mbars) with an interface at 130 mbars and 
14 fixed-p essure levels above (122, 107, 92, 77, 64, 54, 44 , 
34, 24, 14, 7, 3, 1.4, and 0.63 mbars). The horizontal 
resolution of this model is 2’ by 2.5*. The relatively high 
vertical re olution in the upper troposphere and stratosphere 


Table 1. Sources of NO* Used in the Model 

Source Type 

Source Strength 

Reference 

Aircraft 

Fossil fuel burning 
Biomass burning 
Soil microbial activity 
Lightning 
HNO 3 

Transport from the stratosphere 

0.46 

21.0 

6.4 

5.5 

7.0 (2-12.2) 
5.0 
0.4 

iaugheum [19961; Metwally [1995] 
Jenkovitz et al. [1996] 

. iousse et al. [1996]; Atherton [1996] 
'ienger and Levy [1995] 
i Vice et ai [1997a]; Lawrence et al. [1995] 
i ee text 
i se text 


Units are Tg N/yr. 
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together with the advection scheme of Allen et al. [1991] 
assures a numerically accurate representation of tracer 
transport in the lower stratosphere and upper troposphere. In 
particular, the resolution near the tropopause may make this 
one of the better meteorological models available for 
assessing the relative importance of stratospheric sources of 
NO x to the upper troposphere. Convective mass fluxes were 
diagnosed for the same year but were taken from the DAO 
tropospheric model fields (e.g., the GEOS 1 data assimilation 
system; see Schubert et al [1993]), because convective mass 
fluxes were not saved during the assimilation process for the 
stratospheric model. GRANTOUR/CCM 1 , on the other hand, 
uses a highly accurate Lagrangian tracer transport scheme, but 
the basic meteorology driving the model uses a sigma 
coordinate system with 12 vertical levels (mean pressure 
levels 991, 926, 811, 664, 500, 355, 245, 165, 110, 60, 25, 
and 9 mbars) and horizontal resolution of approximately 
4.5’ X7.5\ 

Comparison of model results for 222 Rn from this version of 
the IMPACT model indicates that convection is stronger than 
that in many other models. Zonal average 222 Rn 
concentrations in July reach 10 x 10 ' 21 (vol/vol) in the upper 
troposphere. These predicted upper tropospheric 222 Rn 
concentrations are larger than simulations of 222 Rn that 
directly use the tropospheric GEOS 1 fields [see Allen et al ., 
1996, Allen et al 1997], perhaps because of the differing 
vertical resolution in the tropospheric and stratospheric 
assimilations. The cumulus transport scheme in IMPACT was 
derived from the relaxed Arakawa-Schubert scheme used to 
generate the GEOS 1 meteorological fields. Within the model 
the transport is represented by the cloud transport equation, 

^ = Tooi Ci+l ^* +l ~ qk \~ Ck \^ k ~ qk -' 

where q k is the large-scale tracer mass mixing ratio at level k , 
Q is the tracer mixing ratio within cloud, g is the gravitational 
acceleration (9.8 m s* 2 ), C is the cumulus mass flux (kg m 2 s *) 
derived from the GEOS 1 meteorological fields, and p is the 
pressure in millibars. In these equations, C k +\Q k +[ and C k Q k 
are the transport of tracer into and out of the layer k by cloud 
motions, respectively, while C k +\q k and represent the 

downward transport of tracer by large-scale subsidence 
compensating for the cloud mass fluxes. This equation is 
solved iteratively for each time step together with the 
following equation for Q k , 

(C k +D k )Q k =E k q k +C k+l Q k+ ,, 

where E k and D k are the rates of entrainment and detrainment of 
air into the cloud (kg m 2 s' 1 ), respectively. Below cloud base, 
Q k is set equal to q k , 

Convection in the CCM1 follows a moist adiabatic 
adjustment scheme [Manabe and Holloway , 1975]. The 
convective mass flux was derived from the vertical exchange of 
moist static energy that results from the moist adiabatic 
adjustment, and the following equation for tracer mass was 
solved [Walton et al ., 1988]: 

+ 1 / 2) / Ap k+] )q k+{ - (C(k - 1 / 2)1 Ap k )q k 

dt 100 

where the vertical spatial discretization scheme described by 
Walton et al [1988] has been reversed to use the same 
notation as that described above for the DAO scheme (i.e., 
level *+l is below level k). This equation is solved by using 
an implicit time discretization. The mass fluxes derived from 


the CCM1 can be as much as 20 times higher (on average 2 to 
4 times higher) than those derived in the GEOS-1 assimilation, 
but they are applied to grid-averaged mixing ratios rather than 
to “cloud scale” mixing ratios. Below we show a comparison 
of these two schemes for an idealized experiment. However, in 
spite of the differences in their convective mass transport 
schemes the predicted zonal average upper tropospheric 222 Rn 
concentrations from these two models differ by less than a 
factor of 2. We note that while the convective mass fluxes 
developed for the meteorological fields used by IMPACT are 
more physically based, it is not clear from the radon 
simulations which of the two schemes is the more realistic. 
Thus these two models may represent two plausible extremes 
for vertical and horizontal mixing in the troposphere. 

To carry out a model comparison and separate the effects of 
transport from the effects of chemistry, we have simplified the 
chemistry of NO x and NO,, NO v in the upper troposphere and 
lower stratosphere is the sum of NO x + 2 xN 2 05 + NO 3 +■ 
HO 2 NO 2 + HNO 3 + BrON0 2 + CIONO 2 + PAN plus other minor 
constituents. In the upper troposphere, PAN decomposes very 
slowly and does not play a major role in the cycling of NO x 
versus NO v , so it is not considered further here. Also, in this 
region, BrON0 2 and C10N0 2 are minor components of NO v and 
are also ignored as a source of NO x . The remaining odd- 
nitrogen compounds (N 2 O 5 , N 03 , HONO, HC^NO^ all 
interchange with NO* over the diurnal cycle and are lumped 
into a single tracer family we define as NO,. We define the 
chemistry for NO z from the rates calculated by the Lawrence 
Livermore National Laboratory two-dimensional (2-D) model 
[Kinnison et al . , 1994]. The individual reactions diagnosed 
from the 2-D model that together define the effective loss rates 
for the family NO z are listed below. 

N0 2 + OH -4 HNO 3 
N 2 O s -» 2 XHNO 3 (aerosol) 

C10N0 2 -> HOC1 + HNO 3 (aerosol) 

B 1 ONO 2 -» HOBr + HNO 3 (aerosol) 

N 2 0 5 rainout 
H0 2 N0 2 rainout 

One important process that converts NO- to HNO 3 is 
hydrolysis of N 2 O s occurring on wet aerosols. For the 2-D 
simulation of N 2 0 5 hydrolysis we used an aerosol surface area 
diagnosed from model simulations of sulfate and carbonaceous 
aerosols in the GRANTOUR model [ Penner et al ., 1994; 
Liousse et al 1996]. Because the aerosol reactions mainly 
affect the surface-based sources (which are all continental), we 
used the monthly average aerosol mass concentrations at a 
latitude associated with mainly continental concentrations of 
aerosols, i.e., 22‘E. The model-derived aerosol mass 

concentrations were converted to surface area, assuming a 
lognormal size distribution for ammonium sulfate aerosol that 
had a mode radius for the dry aerosol of 0.05 |im and a 
geometric standard deviation of 2.0 [Kiehl and Briegleb , 
1993]. Carbonaceous aerosol mass from biomass burning, 
fossil fuel burning, and natural organics was converted to 
surface area by using a size distribution and composition 
typical for smoke aerosols [Chuang et al . , 1992]. The dry 
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surface areas for ammonium sulfate aerosol and carbonaceous 
aerosols were increased by a factor of 2.3 and 1 . 3 , 
respectively, to account for aerosol growth to an average size 
at 75% relative humidity. The 2-D rate coefficients, k , were 
then specified according to 

k =— VxSADxv 
4 

where V is the mean molecular velocity, SAD is surface area 
density, and y is the reaction probability, taken as 0 . 1 . This 
expression does not account for gas transport resistance or 
possible room temperature effects [Hu and Abbott, 1997] that 
lower the effective value of y, but sensitivity studies showed 
that decreasing k by as much as a factor of 10 in the lower 
troposphere did not have much effect on calculated NO. in the 
2-D model. Thus the overall rate of this reaction is primarily 
controlled by the rate of the reaction of N0 2 with 0 3 to form 
N 2 0 5 . 

The loss rates developed from the 2-D model in this manner 
were applied in both 3-D models. Figure la shows the 
reciprocal of the total loss frequency (an e-folding lifetime) as 
deduced from the 2-D model. Approximately 30% of the total 
loss rate in the northern hemisphere midlatitude boundary 
layer in July is through the aerosol/N 2 0 5 pathway, decreasing 
to 20% in the upper troposphere. This pathway accounts for 
60% of the loss in the winter hemisphere. We also used the 2-D 
model to diagnose the ratio of NO* to NO, (shown in Figure lb) 
for each month. This ratio was used to diagnose NO r given the 
predicted zonal average concentrations of NO z from each 3-D 
model. 

This procedure restricts our calculated NO* to that obtained 
from the primary sources without recycling from HN0 3 The 
secondary source of NO* from HN0 3 recycling was included 
separately by diagnosing the loss frequency for HN0 3 
photolysis and its reaction with OH to form N0 3 in the 2-D 
model. This approach allows us to examine the range in 
background tropospheric NO* concentrations resulting from 
uncertainties in concentrations and recycling rates of HNO 3 . 
This scheme differs from that taken in the recent study by 
Lamarque et aL [1996], who tagged each nitrogen-containing 
molecule from each source of NO* in their model and allowed 


feedback between the calculated NO*, other nitrogen 
reservoirs and the odd hydrogen system. Our approach 
overestim ites the increase in NO* that results from a given 
increase n source strengths because we have neglected 
feedbacks between the odd hydrogen abundance and the 
concentration of NO*. Figure 2 , for example, shows how the 
lifetime cf NO* varies with NO* concentration at different 
altitudes i 1 the midlatitude troposphere. While the dependence 
is weak in winter, in summer, below about 500 pptv NO*, the 
concentrat on of OH (which mainly determines the lifetime of 
NO*) is primarily controlled by the reaction of HO 2 with NO. 
Thus, in this concentration regime, increases in NO* sources 
and concentrations lead to increases in OH, thereby decreasing 
the lifetime of NO*. Above the 500 pptv level, further 
increases n NO* increase the lifetime of NO* because the 
reaction of OH with N0 2 effectively removes odd hydrogen so 
that OH decreases as NO* increases. Our linearized chemical 
scheme will overestimate the sensitivity of NO* concentration 
to changes in source strength when NO* is less than 500 pptv 
and when recycling from HN0 3 is calculated separately. On the 
other hand, this simplification allows us to separately explore 
the effect of recycling and also has the advantage of allowing a 
direct comparison of the effects of different representations of 
transport on the NO* concentration without the complicating 
influence of changes in the NO* lifetime. Furthermore, as we 
show below, model differences in transport can lead to far 
larger differences than the factor of 2 change in lifetime shown 
in Figure lb. A similar linear approach was also used in the 
study by K^aus et al. [1996] and Kohler et aL [1997]. 

Table 1 describes the NO* sources used in this study together 
with literature estimates for the range of source strengths 
considered here for the lightning source. The table also shows 
the derived source strength from a background concentration of 
100 pptv HNO 3 . The surface sources noted in the table were 
input into :he lowest model layer for the IMPACT model and 
into the 1(00-900 mbar domain of GRANTOUR. The aircraft 
emissions 1 sed in these calculations are based on air traffic for 
each montf of 1992. Emission inventories of NO* for both 
scheduled [Baughcum et al ., 1996] and nonscheduled 

[Metwally, 1995] air traffic were combined on a 1° X 1 ° X 1 
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Figure 2. Lifetime of NO, as a function of NO, concentration in (a) winter northern midlatitudes and (b) 
summer midlatitudes. 


km (pressure altitude) grid and then interpolated onto the 
model grids. The lightning sources were developed from 
monthly mean NO* production rates per unit surface area as a 
function of geographic region based on cloud top heights 
derived by the International Satellite Cloud Climatology 
Project (ISCCP) and defined at three levels: 115, 245, and 375 
mbar [Price et al , 1997a). Below cloud top, lightning NO, was 
distributed with altitude according to an approximate C-shaped 


profile developed for deep convective events for the tropics, 
midlatitudes, and oceans (Table 2, from Pickering et al. [1996] 
and K.E. Pickering et al. (Vertical distributions of lightning 
NO, for use in regional and global chemical transport models, 
submitted to Journal of Geophysical Research, 1998). Profiles 
in the table were scaled to the ISCCP cloud top heights. 
Originally, lightning NO, was further restricted to levels 
below the tropopause defined from the temperature profile in 


Table 2. Percentage of Mass of NO, Produced by Lightning as a Function of Altitude for Deep Convective 
Clouds in Three Regimes — 


Height, km 


Tropical Marine* Tropical Continental Midlatitude Continental 


0-1 

0.7 

1-2 

0.5 

2-3 

0.7 

3-4 

0.7 

4-5 

0.6 

5-6 

0.7 

6-7 

1.3 

7-8 

4.8 

8-9 

12.1 

9-10 

19.2 

10-11 

20.7 

11-12 

19.7 

12-13 

11.5 

13-14 

4.1 

14-15 

1.3 

15-16 

0.5 

16-17 

0.3 

17-18 

0.3 

18-19 

0.2 

19-20 

0.2 


3.4 

21.9 

1.4 

4.4 

1.5 

1.2 

1.6 

18 

1.1 

3.5 

1.3 

4.4 

1.6 

3.9 

2.8 

4.7 

4.6 

5.9 

6.5 

7.8 

8.4 

10.3 

11.0 

10.8 

14.8 

9.3 

17.2 

6.6 

13.6 

2.7 

6.4 

0.6 

2.1 

0.1 

0.6 

0.0 

0.1 

0.0 

0.0 

0.0 


*A11 marine sources of lightning NO r in the model used the “tropical marine" profile developed by Pickering et al 
[1996] and K.E. Pickering et al., (Vertical distributions of lightning NO t for use in regional and global chemical transport 
models, submitted to Journal of Geophysical Research , 1998). 
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each model. However, this procedure sometimes injected NO v 
into the a region with largely stratospheric air, and we 
therefore restricted cloud top lightning NO* to be below the 
tracer tropopause defined in section 3. 

3. Stratospheric Source of NO* 

We expect that the differentiation between air with 
substantial concentrations of stratospheric tracers and air with 
mainly '’tropospheric composition" should be represented in 
the model by a "tracer tropopause." Air parcels above this 
surface should have longer exchange times with the lower 
troposphere, while parcels below this level communicate more 
rapidly with surface and boundary layer air. We have 
investigated the location of this surface in both GRANTOUR 
and IMPACT by an idealized experiment in which we injected 
NO_ v into the stratosphere uniformly between 50 mbar and 1 0 
mbar at the rate of 0.4 Tg N/yr and removed it below 800 mbar 
with a 5-day time constant (case 1) and with a 1-hour time 
constant (case 2). Regions in which the case 1 and 2 tracers 
differ by more than 5% clearly define regions that are 
relatively strongly coupled to the lower troposphere. We 
define the surface between this model-defined region of 
tropospheric air and air with substantial stratospheric 
influence as the "tracer tropopause." This level is 

approximated in midlatitudes in each model by the 100 pptv 
stratospheric tracer concentration surface. The use of the 1 00 
pptv surface as a boundary allows us to cleanly define air that 
has substantial stratospheric influence within the two models. 
The location of this surface at each grid point is shown in 
Figure 3 for January 1 for the two models. As Figure 3 shows, 
this is a highly dynamic surface that changes with longitude as 
well as time (the vertical range of points with 100 pptv 
concentration increases by about a factor of 2 if a whole month 
of data is plotted). The heavy line in this figure and in 
subsequent figures shows the monthly average model- 
calculated tracer tropopause in January and July. In the tropics 
this surface is nearly identical in the two models. North of 
40°N in July and north of 30°N in January, however, the 
GRANTOUR model’s tropopause begins to descend more 


rapidly th in that in IMPACT. At the north pole, GRANTOUR 
has a tracer tropopause level of approximately 400 mbar, 
while IMPACT has a tracer tropopause level of 300 mbar. It is 
clear that the effective lifetime for NO v from aircraft sources 
that are largest between 200 and 300 mbar and between 40°N 
and 60°N will differ significantly between the two models as a 
result of he placement of this source in relation to each 
model’s tracer tropopause. In addition, the contribution of 
lightning io NO, could be significantly larger if the source is 
placed abt ve this level. Similar model differences occur in the 
southern hemisphere poleward of 30° latitude. 

We no e that the choice of the 100 pptv level of 
stratosphe ic NO v tracer as a surrogate for the 5% difference 
level between our two tracer experiments is significantly 
smaller than the concentration of NO v normally associated 
with the physical tropopause (defined by the level above 500 
mbar at which the lapse rate becomes smaller than 2°K km' 1 ), 
and in fact, our 100 pptv level lies below each model's average 
physical tropopause in midlatitudes. We note that the physical 
tropopause shows similar scatter about the average, similar to 
our 100 pptv level. According to observations, the average 
value of NO v near the physical tropopause is close to the 500 
pptv concentration level but may vary between 200 and 800 
pptv [Murphy et al 1993]. Data presented by Murphy el al. 
[1993] also show a decrease to an average value of 250 pptv 
within 50G m of the physical tropopause (see Figure 10 of that 
paper). In both models we find that the zonal mean 250 pptv 
NO v contour in midlatitudes approximates the zonal mean 
physical tropopause, while the 500 pptv level is above the 
average physical tropopause in midlatitudes. 

A number of reasons may explain why our model-derived 
tracer tropopause NO v value does not agree with the 500 pptv 
level of Murphy et al. [1993], First, we are measuring the 
surface that describes reasonably fast (5 to 10 days) 
communication with air below 800 mbar, a surface that 
describes where sources may build up in the model to 
concentrati >ns higher than those they would attain within the 
lower trope sphere. Second, even our high-resolution model, 
IMPACT, 1 as a vertical resolution in this region of only about 
1500 m. Thus it is perhaps remarkable that the models are able 



Figure 3. Instantaneous vertical location of the 100 pptv contour >n January 1 for the stratospheric tracer 
experiment defined in section 3 for (a) GRANTOUR and (b) IMPACT. 
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to come within a factor of 2 of the value of NO v at the physical 
tropopause, given that real atmospheric concentrations may 
change by that much in only 500 m. Thus numerical diffusion 
may be causing the concentrations in our modeled tracer- 
tropopause to be less than they would be in a higher-resolution 
model. Third, if one examines plots of NOy, 0 3 , or their ratio 
when plotted against pressure level (as we do approximately, 
by plotting NO* results in sigma and eta coordinates), the 
features seen in the plots may not precisely correspond to 
features seen in data when plotted as distance from the physical 
tropopause. For example, Figure 6 of Murphy el ai [1993] 
shows that N0*/0 3 is approximately constant above about 100 
mbar in midlatitudes and is equal to 0.003, but N(y0 3 
increases to 0.007 at 300 mbar, a level consistent with an 
average of 100 ppbv 0 3 (the value usually associated with the 
tropopause). Thus the level that one might associate with the 
tropopause on the basis of the 0 3 concentration has a NQ/O3 
ratio outside the range usually considered stratospheric. 
Clearly, local variability in the pressure level of the physical 
tropopause can make analysis of precise concentrations of 
tracers difficult in this coordinate. Finally, it is possible that 
lightning-produced NO v in the upper tropical troposphere 
could enter the tropical stratosphere, supplying the lower 
stratosphere with more NO v than one would calculate in a 
model experiment that includes only the stratospheric source. 
This NO v would be a minor component of total NO* at higher 
altitudes in the stratosphere but may contribute to the 
concentration of NO v at the tropopause. 

The model study undertaken here must differentiate between 
primary NO* as that which is directly input from a given source 
and secondary NO* as that which is recycled from HN0 3 (other 
sources such as PAN, as noted above, are not important in the 
upper troposphere and are not included). This problem is not 
easily formulated with respect to the stratospheric source of 
primary NO*, since most NO v in the stratosphere is HN0 3 and 
less than 10% enters the troposphere as primary NO*. We 
present a new approach here that clearly separates primary and 
secondary NO* from the model stratosphere. As we noted 
above, the "tracer tropopause" is associated with air having 
greater than 100 pptv of stratospheric NO v in both models. 
Therefore, in our simulations of the stratospheric source of 
primary NO*, we defined all air that had a stratospheric NO v > 
100 pptv as "stratospheric." Within this region of the model, 
both NO* and HN0 3 were carried as prognostic variables, so 
that NO* formed HN0 3 with a time constant defined by the 2-D 
model (i.e., by the reactions that form HN0 3 ) and HN03 was 
recycled back to NO* with a time constant defined by the 
photolysis of HNO3 and its reaction with OH as diagnosed 
from the 2-D model. When NO v (= NO* + HN0 3 ) was below 1 00 
pptv, we followed the procedure outlined in section 2 for the 
tropospheric sources of NO*, namely, NO* was removed at the 
local rate determined from the 2-D model (Figure la) and there 
were no other sources of NO* (no recycling from HN0 3 to NO*, 
which is treated separately; see section 4). 

Figure 4 shows the zonal average concentrations of primary 
stratospheric NO* from the GRANTOUR and IMPACT models. 
This figure demonstrates that mixing of this stratospheric 
tracer below the tracer tropopause is, at least at northern 
midlatitudes, surprisingly similar in the two models. A strong 
downward mixing near 30°N is evident in the IMPACT model, 
especially in July. We note that the stratospheric NO* mixing 
ratios appear quite flat in both models throughout the tropical 
upper troposphere, a region normally associated with upward 


movement of air. The GRANTOUR mixing ratios in January 
even appear consistent with a downward moving source in this 
region. The flat contours in the tropical upper troposphere are 
associated with horizontal mixing in both models. To 
investigate the source of the apparent downward movement of 
tracer in the January simulations in GRANTOUR, we performed 
a model experiment in which all mixing processes within 
GRANTOUR were turned off. Thus the only transport was that 
due to large-scale advection. The appearance of the source near 
the tropical tropopause was diminished because downward and 
horizontal transport by mixing was zeroed, but the feature was 
still present at lower altitudes. We compared these features 
with the zonal average vertical velocities from the model. This 
comparison showed that the apparent downward moving source 
in the tropics is actually caused by downward motions in the 
subtropics. This is followed by horizontal advection, which 
makes it appear (in a zonal average plot) as though the tracer is 
originating from the stratosphere in downward motion at the 
tropics, while in fact this downward motion occurs at different 
latitudes. Similar processes take place in the July simulations 
but do not exhibit themselves as a possible downward moving 
stratospheric source in the tropics. 

The direct injection of NO* into the midlatitude troposphere 
below the tracer tropopause from the stratosphere is quite 
small in both models. Typically, the free tropospheric values, 
well away from the tracer tropopause, are less than 5 pptv. The 
most abundant form of NO v from stratospheric injection, 
however, is predicted to be HN0 3 , with typical concentrations 
of 20 to 50 pptv in the free troposphere in both models (not 
shown). HN0 3 concentrations as large as 100 pptv could only 
produce an additional 5-10 pptv of NO* in the troposphere, as 
is shown along with the other NO* sources in section 4. 
Therefore the total NO* concentration due to stratospheric 
sources (i.e., derived from a combination of direct NO* 
transport and HN0 3 chemical conversion to NO*) is at most of 
the order of 20 pptv below the tracer tropopause. Above the 
tracer tropopause the total NO* was calculated from the 
exchange with HNO3. The values calculated in July north of 
45 # N and between 200 and 300 mbar are between 5 and 50 pptv 
in GRANTOUR and between 2 and 50 pptv in IMPACT. They 
are even smaller in the winter upper troposphere and lower 
stratosphere. These concentrations are considerably smaller 
than the measured concentrations of 120-150 pptv NO and 
100-300 pptv NO made off the coast of North America during 
the June 1984 Stratospheric Ozone Experiment (STRATOZ III) 
[Drummond et ai . , 1988] and during the January 1991 
Tropospheric Ozone Experiment (TROPOZ II) [ Rohrer et al . , 
1997], respectively. They are also far smaller than most the 
upper tropospheric measurements of NO and NO* compiled by 
Emmons et al. [1997]. This finding argues against these 
measured concentrations being explained, on average, by a 
stratospheric source. 

4. Contribution of Tropospheric Sources to Upper 
Tropospheric NO* 

Figure 5 shows the zonal average concentrations of NO* for 
the six tropospheric sources treated here from the IMPACT 
model in July, while Figure 6 shows those from GRANTOUR. 
We have superimposed on these graphs the average tracer 
tropopause from Figure 3. The altitude and magnitude of peak 
concentrations for both models from the aircraft source 
(Figures 5a and 6a) are quite similar. Note, however, that the 
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Figure 4. Predicted zonal average primary NO x mixing ratio from the stratospheric source (pptv) in (a), (c) 

January and (b), (d) July in the GRANTOUR and IMPACT models, respectively. Superimposed on these plots is 
the tracer tropopause height in each of the models. The troposphere is defined as the region where the 
stratospheric tracer concentration was less than 100 pptv. 

aircraft source is deposited mainly within the model differences in the upper troposphere, the surface 

stratosphere for GRANTOUR, while it is partly in the concentrati ms are similar in the two models. This large 

stratosphere and partly in the troposphere for IMPACT. This difference h the contribution of surface sources to the upper 

difference, together with the more rapid vertical mixing in the troposphert was unexpected, since the model -predicted radon 

IMPACT model, brings NO* from the upper level source down concentrati ms were within a factor of 2. However, NO as 
to the surface more rapidly than happens in GRANTOUR. defined hers has a much shorter lifetime at the surface in 

Consequently, the 5 and 2 pptv concentration contour comparison with radon. The ^-folding lifetime for radon (5.5 

intervals from aircraft NO r reach farther into the troposphere days) allows venting of boundary layer air by convective 

in IMPACT in comparison with the same contour intervals in mixing processes to compete effectively with its chemical 

GRANTOUR. Despite this rapid mixing the IMPACT model has removal in mth models. Here the e-folding lifetime for NO. is 

a peak zonal average NO* from aircraft that is >20 pptv, a only 1 day near the surface, effectively shutting off vertical 

result of the higher resolution in this model. transport ftr the slower overturning in the GRANTOUR model 

Figures 5b and 6b, 5c and 6c, and 5d and 6d depict the zonal but allowirg some emissions to be carried to the upper 

mean NO x concentrations from biomass burning, fossil fuel troposphere in the IMPACT model. 

burning, and the soil source, respectively. The upper Concenti ations of NO x from lightning (source rate of 7 Tg 

tropospheric NO x concentration from these surface sources are N/yr) are shown in Figures 5e and 6e. In the upper troposphere, 

very different in the two models. Thus concentrations in both models predict concentrations reaching 100 pptv. In 

IMPACT reach 20 pptv, 100 pptv, and 50 pptv for the three preliminary simulations we had inadvertently placed the 

surface sources in the upper troposphere in July, respectively, lightning source in the IMPACT model above the tracer 

while the corresponding concentrations from GRANTOUR are tropopause, causing a factor of 2 higher peak zonal average 

only <1 pptv, 5 pptv, and 2 pptv. In contrast to these concentratio is of NO x [Penner et al., 1996)]. Here we see that 
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Figure 5. Predicted zonal average NO* mixing ratio (pptv) in July from the IMPACT model from (a) aircraft, 
(b) biomass burning, (c) fossil fuel burning, (d) soils, (e) lightning, and (f) the HNO 3 source. 


while peak concentrations are similar, just as for the aircraft 
source, the stronger vertical mixing in IMPACT brings more 
lightning NO* toward the surface. 

Finally, Figures 5f and 6 f show the NO* concentration that 
results from the specified background source from HNO 3 , 


assuming a concentration of 100 pptv HNO 3 . Again, the most 
important difference between the two models appears to be the 
rate of vertical mixing, with NO* concentrations from HNO 3 
penetrating farther toward the surface in IMPACT than in 
GRANTOUR. It can also be noted that this is a relatively minor 
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Figure 6. Predicted zonal average NO, mixing ratio (pptv) in Ju y from the GRANTOUR mode! from (a) 
aircraft, (b) biomass burning, (c) fossil fuel burning, (d) soils, (e) ligitning, and (0 the HN0 3 source. 


source of NO* in the upper troposphere in comparison with 
that from lightning in both models or that from fossil fuels in 
IMPACT. As we noted previously, this background source 
contributes 5-10 pptv NO* in the middle latitude upper 
troposphere and only reaches 20 pptv in the tropical upper 


troposphere. Thus it is unlikely that upper tropospheric NO* is 
largely fron the stratosphere. 

We inv< stigated the effect of the convective mixing 
algorithms in the two models with an off-line single-column 
experiment In this experiment an initial profile for an 
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idealized tracer was specified wherein the initial mixing ratio 
was 1 below a sigma level of 0.74 and zero above that level. 
We evaluated the mixing ratio after a 48-hour simulation in 
which the tracer was redistributed according to the convective 
mixing algorithm in each model. We examined two cases: one 
in which the lifetime of the tracer was 1 day below the 0.74 
sigma level and one in which the lifetime was 5 days below 
0.74 sigma. In both cases there was no loss above 0.74 sigma. 
An idealized convective mass flux profile was assumed in units 
of kg m' 2 s' 1 , which had values of simply the sigma coordinate 
divided by 100. So, for example, the convective mass flux at 
0.58 sigma was 0.0058 kg m' 2 s' 1 . The convective mass fluxes 
between sigma levels of 0.2 to 0.1 were then forced to decrease 
linearly to zero. There was no mass flux for sigma levels 
above 0.1. An exact comparison with the model results shown 
in Figures 5 and 6 is problematic, of course, because the 
magnitude of the derived convective mass fluxes from the two 
models differ by up to a factor of 20, but these idealized 
experiments serve to illustrate what the effect of the two 
algorithms is. To approximate the effect of different mass 
fluxes, we applied the algorithm used in GRANTOURyCCM 1 , 
using the same values of convective mass flux, as well as 2 
times and 5 times these values. Figure 7 shows the result. The 
tracer profile developed by using the IMPACT algorithm 
shows apeak in mixing ratio above the 0.25 sigma level in 
each case, while the GRANTOUR profiles decrease smoothly 
for the 5-day tracer lifetime and develop a peak in 
concentration at lower altitudes in the 1-day tracer lifetime 
experiment. For larger convective mass flux values the mixing 
ratios for GRANTOUR increase but never get as high as those 
predicted with the IMPACT algorithm at the highest altitudes. 
The redistribution of tracer mass by convection clearly 
depends also on the initial profile used in this type of 
experiment. A similar off-line calculation in which the tracer 


was initially well mixed vertically showed that the 
GRANTOUR algorithm with 5 times the convective mass flux 
places more material above the 0.25 sigma level than does the 
IMPACT algorithm. Figure 7 shows that the schemes in both 
IMPACT and GRANTOUR can take material from the boundary 
layer into the free troposphere equally well on a 1-day 
timescale. The amount actually delivered to the upper 
troposphere depends on the statistics of convective processes 
in the models, the magnitude of the mass flux, and the profiles 
that develop within the model. Figure 7a, however, 
demonstrates the fact that the treatment of convection can lead 
to large differences in upper tropospheric NO* concentrations 
if the lifetime of the tracer is only of the order of 1 day. This 
comparison demonstrates the need to evaluate the model 
treatment of convection, using a species that is much shorter- 
lived than the conventional 222 Rn to evaluate the model 1 s 
capability to simulate NO*. 

Figures 8 and 9 show the same results as Figures 5 and 6 , 
respectively, but for the month of January. Concentrations 
from aircraft are somewhat lower in January than in July. On 
the basis of chemical removal rates, which are more rapid in 
summer, we would expect higher concentrations in January. 
Thus the aircraft NO* concentrations are most likely controlled 
by seasonal differences in the dispersion by the 
meteorological fields together with the seasonal variation in 
aircraft emissions, which increase by about 30% between 200 
and 300 mbar and 40° to 60’ N from January to July. 

Figures 8b and 9b may be compared with Figures 5b and 6b, 
respectively, to compare the biomass burning sources in 
January and July. As these sources are mainly tropical, and 
both photochemistry and vertical mixing are similar in both 
seasons in the tropics, the distributions are similar, with July 
peaks occurring south of the equator and January peaks 
occurring north of the equator, consistent with the biomass 



Figure 7. Predicted mixing ratio for an idealized experiment using two convective mixing schemes. 
IMPACT results are shown by triangles. Squares show the results from the algorithm in GRANTOUR using the 
same convective mass flux as IMPACT, while diamonds and circles used 2 times and 5 times the mass flux, 
respectively. The lifetime of the tracer was (a) 1 day and (b) 5 days below the 0.74 sigma level. 
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Figure 8. Predicted zonal average NO, mixing ratio (pptv) in Jan rnry from the IMPACT model from (a) 
aircraft, (b) biomass burning, (c) fossil fuel burning, (d) soils, (e) ligh ning, and (f) the HN0 3 source 


source distributions. The distribution of fossil fuel NO*, on July at no them midlatitudes at the surface the strong 

the other hand, is quite different between the two seasons convection in IMPACT brings peak concentrations of over 

(compare Figures 8c and 9c with Figures 5c and 6c, 100 pptv io the upper troposphere (Figure 5c), while in 

respectively). Despite the stronger photochemical removal in January (Fij ure 8c) the concentrations are a factor of 10 lower 
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Figure 9. Predicted zonal average NO, mixing ratio (pptv) in January from the GRANTOUR model from (a) 
aircraft, (b) biomass burning, (c) fossil fuel burning, (d) soils, (e) lightning, and (0 the HNO 3 source. 


at the same altitudes. The slower vertical mixing in the 
GRANTOUR model does not introduce a large seasonal 
variation in the resulting upper tropospheric NO,. 

The predicted concentrations from the soil source at 
northern midlatitudes are higher in July both at the surface and 


at upper tropospheric levels, reflecting emission rates that are 
higher in the midlatitude northern hemisphere in summer. The 
HNO3 source for January is a near reflection (about the equator) 
of the source for July, although there appears to be more 
vertical mixing in January in the southern hemisphere in the 
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IMPACT model than there is in July in the northern constnicle 1 for the IMPACT model are considerably higher i n 

hemisphere. The opposite appears to be the case for the the upper troposphere/lower stratosphere than are the upper 

GRANTOUR model (less mixing in the southern hemisphere and lower imits constructed for GRANTOUR. As we noted from 

summer than in the northern hemisphere summer). The our previo js examination of different source strengths, this 

distribution of NO, from lightning is similar in the two result is m tinly due to the stronger vertical mixing present in 

seasons except that the stronger source in July in Northern the IMPACT model and its transport of surface sources. The 

midlatitudes is evident in the large area surrounded by the 50 upper and lower limits for NO, vary by about a factor of 5 for 

Pptv contour in the IMPACT model and the 20 pptv contour in the GRAN TOUR model, while the range is only about a factor 

the GRANTOUR model. of 2 in IMPACT. In keeping with this result the contribution 

to upper tropospheric NO, from all anthropogenic sources is 
5. Conclusion much ,ar 8 < r in IMPACT (ranging from 40% to 60% in northern 

midlatitudi s in July) than it is in GRANTOUR (10% to 30%). 

The upper and lower limits of total NO, (all seven sources) These emulations demonstrate the large uncertainties 
in January and July are shown in Figures 10 (IMPACT) and 1 1 associated with prediction of NO, in the upper 

(GRANTOUR). The upper limit assumes 12 Tg N/yr of troposphei e/lower stratosphere. Of the three main 

lightning NO, and a background of 200 pptv HNOj, whereas determinants of upper tropospheric NO, (transport, lightning 

the lower limits assume 2 Tg N/yr of lightning NO, and 50 source strength, and recycling from HNO3) that hamper 

pptv HN0 3 In these calculations the source of NO, from HNO3 abilities to predict upper tropospheric/lower stratospheric 

was added below the model's tracer tropopause, while above NO,, the 1 ghtning source strength and vertical transport are of 

that level, NO, from the stratospheric tracer experiment most importance. Predictions of midlatitude upper 

(Figure 4) was used. Both the upper and lower limits tropospheric NO, of the order of 100 pptv in July are 



Figure 10. Predicted zonal average NO, mixing ratio (pptv) in IMPACT in (a), (b) January and (c), (d) July 
using the upper and lower bounds for lightning (12 and 2 Tg N/yr) and HNO3 (200 and 50 pptv background 
concentrations) sources, respectively. 
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Figure 11. Predicted zonal average NO* mixing ratio (pptv) in GRANTOUR in (a), (b) January and (c), (d) 
July using the upper and lower bounds for lightning (12 and 2 Tg N/yr) and HNO3 (200 and 50 pptv 
background concentrations) sources, respectively. 


consistent with either slow vertical mixing and a large 
lightning source or more rapid vertical mixing and a smaller 
lightning source (compare Figures 11c and lOd). Thus, to 
improve model assessment of ozone change in the upper 
troposphere/lower stratosphere, it is most important to 
improve the vertical transport treatment in current models. 
Further, to test model treatment of convection and vertical 
transport as it applies to the NO* problem, it is important 
either to find a tracer whose lifetime is more similar to that of 
NO* or to use methods other than comparison of monthly 
average model fields with data for testing the models (see 
below). While this approach is useful for many longer-lived 
species, comparison of predicted monthly average 222 Rn 
concentrations with data does not provide a sufficiently 
accurate test for the short-lived NO*. Comparison of average 
H2O and cloud amounts with data (as used by the climate 
modeling community) would also be insufficiently accurate for 
testing the transport of short-lived NO*. 

Evaluation of the effects of aircraft emissions present a 
particularly difficult model challenge. To understand the model 
response to these emissions, it is very important to understand 
the model-derived tracer tropopause. Some models may place 


aircraft emissions entirely within regions with large influence 
of stratospheric air, while others place them in regions that 
exchange relatively rapidly with air that has only a 1-day 
lifetime for NO*. Entirely different concentration profiles (and 
perturbations to ambient NO* and other species) are expected 
as a result of these differences. 

Future verification of any CTM’s prediction of tropospheric 
NO* and hence O3 perturbations due to increases in 
anthropogenic NO* (e.g., aircraft, surface combustion) requires 
that each component of NO* be accurately modeled as well as 
the total. At present, measurements give us only the total NO*, 
and we must rely on model intercomparisons and independent 
evaluations of individual source strengths to test different 
CTMs. One method that seems promising, however, is to use 
NO*/NO v as a diagnostic for fresh emissions and correlations 
of different species with NO* or NO v to diagnose the frequency 
with which air parcels with different sources of NO* contribute 
to upper tropospheric air. For example, correlation of NO v 
with O3 and anticorrelation with H2O can be used to diagnose 
stratospheric air, and correlations of NO* with NMHCs, CO or 
222 Rn can diagnose air that had recently come from surface 
sources. These tracer/tracer correlations together with 
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characterization of the ratios of the change in NO v with change 
in CO or C0 2 have been used to identify air perturbed by recent 
aircraft emissions as well as other sources [Zheng et al., 
1996]. Testing the ability of models to reproduce such 
correlation frequencies appears to be one of the more 
promising methods for gaining confidence in such models and 
for better quantifying the lightning source. 
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